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ABSTRACT 


Two  new  ballistic-type  self-dispersing  bomblets  which  utilize  passive 
control  of  the  radial  orientation  of  the  trimmed  lift  force  are  described 
and  the  motion  of  these  bomblets  as  determined  by  simplified  theory  and  bv 
exact  six-degrees-of-freedom  simulation  is  discussed.  Some  of  the  unique* 
aerodynamic  characteristics  of  these  bomblets  as  determined  from  wind  tunnel 
tests  .'•.re  also  described.  Impact  pattern  predictions,  using  the  Monte  Carlo 
method,  show  that  properly  designed  S-Curve  and  roll -through -zero  bomblets 
can  achieve  large  uniform  patterns  for  low-altitude  high-speed  delivery 
conditions;  however,  configurational  asymmetries,  if  not  suitably  restricted, 
can  significantly  modify  the  bomblet  performance. 
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SECTION  I 


INTRODUCTION 


The  inr.reasing  utilization  of  cluster -type  weapons  for  tactical  war¬ 
fare  has  led  to  the  need  to  consider  various  bomblet  dispersion  techniques. 
During  the  past  several  years  two  new  self-dispersing  bomblet  concepts  which 
utilize  the  lift  force  resulting  from  body  incidence  as  a  means  of  obtain¬ 
ing  dispersion  have  been  investigated.  Each  employs  a  novel  means  of 
passively  restricting  the  rotation  of  the  lift  vector  about  the  flight  path 
as  a  means  of  greatly  increasing  the  dispersion  potential  in  the  presence 
of  configurational  asymmetries. 


The  S-Curve  bomblet  is  comprised  of  an  axially  symmetric  body-fin 
configuration  which  is  designed  to  provide  an  unstable  restoring  moment  at 
small  angle  of  attack  and  a  stable  pitching  moment  slope  at  a  large  trim 
angle  of  attack.  The  S-Curve  name  is  derived  from  the  shape  of  the  non¬ 
linear  pitching  moment  curve.  The  salient  feature  of  the  S-Curire-type 
bomblet  is  that  the  trim  incidence  is  independent  of  body-roll  attitude 
such  that  the  direction  of  the  lift  force  vector  with  respect  to  a  vertical 
plane  is  either  fixed  or  slowly  changing,  even  though  the  bomblet  has  an 
accidental  spin.  The  use  of  this  idea  for  bomblet  dispersion  was  first 
described  by  Merchant  and  Pope.(l) 

The  roll -through- zero  bomblet  is  comprised  of  a  lifting-body  configu¬ 
ration  equipped  with  a  roll  producing  device  which  will  cause  the  bomblet 
to  roll  in  a  direction  opposite  to  its  initial  direction  of  roll.  The 
dispersion  of  this  type  bomblet  depends  upon  the  temporary  and  limited  change 
in  the  orientation  of  the  lift  vector  as  the  bomblet  passes  through  zero 
roll.  The  fact  that  large  dispersion  can  result  from  a  roll -through-zero 
motion  has  been  noted  earlier  with  respect  to  the  development  of  re-entry 
vehicles. (2^ 

This  report  describes  the  analytical  dispersion  predictions  for  several 
candidate  S-Curve  and  roll -through- zero  bomblet  configurations. 
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SECTION  II 


FLIGHT  DYNAMICS 


S-Curve  Bomb let 


Aerodynamic  Characteristics.  Two  considerations  dominate  the  design 
of  the  S-Curve  bomblet:  first,  the  attainment  of  a  proper  nonlinear  restoring 
moment  characteristic  throughout  the  operating  Mach  number  range;  second, 
restriction  of  roll  dependent  aerodynamic  forces  and  moments,  particularly 
the  side  moment,  which  causes  undesirable  coning  motion.  The  latter  is 
accomplished  by  increasing  the  axial  symmetry  (e.g.,  number  of  fins)  and 
by  limiting  the  trim  angle  of  attack.  The  S-Curve  moment  characteristics 
are  obtained  most  effectually  by  a  combination  of  nose  rounding,  boattailing, 
and  utilization  of  very  low  aspect-ratio  fins  or  sub-caliber  stabilizers. 

Figure  1  illustrates  two  basic  S-Curve  bomblet  designs  which  have  been 
extensively  investigated  in  the  wind  tunnel. Both  are  four  calibers 
in  length;  however,  the  first  has  a  cylindrical  afterbody  while  the  second 
has  a  boattail  afterbody.  Figure  2  illustrates  the  pitching  moment,  normal 
force,  and  trim  characteristics  of  the  two  basic  S-Curve  configurations. 

The  two  S-Curve  bomblet  configurations  differ  essentially  in  their  transonic 
trim  characteristics.  At  transonic  Mach  number  the  cylindrical  afterbody 
bomblet  exhibits  large  variations  in  trim  angle,  while  for  the  bomblet  with 
the  boattail  afterbody  the  trim  angle  is  nearly  invariant  with  Mach  number. 

The  induced  roll  torque  and  side  m.oment  coefficients  at  $  =  15  degrees 
(approximately  maximum  values),  as  well  as  the  magnus  moment  coefficient, 
are  plotted  in  Figure  3  as  a  function  of  angle  of  attack  for  the  cylindrical 
afterbody  configuration.  Two  different  magnus  moment  curves  are  plotted; 
the  original  estimate  is  based  on  test  results  of  an  Australian  Defence  ^.g.^ 
Scientific  Service  Weapons  Research  Establishment  designed  S-curve  bomblet,'  ' 
and  is  used  for  some  flight  dynamic  calculations  which  will  be  described 
later.  All  of  the  moments  are  high  nonlinear  functions  of  angle  of  attack; 
however,  for  trim  angles  of  attack  between  8  and  16  degrees  the  moments,  with 
the  exception  of  the  side  moment,  are  small. 


Ca)  S-Curve  Bomblet:  Cylindrical  Afterbody 


(b)  S-Curve  Bomblet:  Boattail  Afterbody 


Cc)  Roll -Through-Zero  Bomblet 


Figure  1.  Bomblet  Configurations 
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Zero  Coning  Theory.  The  tendency  for  the  angle  of  attack  plane  to 
remain  stationary  (i.e.,  for  the  coning  to  be  zero)  is  the  fundamental 
reason  for  the  large  dispersion  of  the  S-Curve  bomblet.  The  requirements 
for  zero  coning  can  be  stated  precisely  by  considering  the  differential 
equation  for  the  variable,  (^,  which  is  the  orientation  of  the  angle  of 
attack  plane  with  respect  to  a  vertical  reference  plane. 


-  Pi  g  (Xj.  sin  Y  I  Cjj  p  S  V  n  p  S  d^  V  otj."!  ^ 

a  _,  ({»  -  -  -  - +  -  <j) 

L  V  2m  J 

g  p  S  d^  cos  Y  1  ~ 


P  8  cos  yI  ~  [  g  p  S  d^  cos  Y 

+  -  jeos  4i  -  I  -  sin  (i»  = 

L  V  J  L  4  J 


(1) 


p  S  d  V 


CSM  0  S  d 


2m 


The  existence  of  a  stj|ady-state  solution  for  in  Equation  (1)  is 
determined  by  setting  ^  =  0.  This  results  in  Equation  (2), 


I  7  C„  p  S  V 

I  g  cos  Y  cos  <t>  _  N 


2m 


g  (Xj.  sin  Y 


h^  g  p  S  d  cos  Y  sin  ^ 


.  1  Si  p  p  s  V 
4  p 


(2) 


-  ^  P  V’  S  d  =.  0 

where  the  first  term  Cin  braces)  is  the  gyroscopic  moment  due  to  rate  of 
change  of  flight  path  curvature,  the  second  term  is  the  damping  moment 
due  to  circular  motion,  and  the  third  and  fourth  terms  are  the  magnus  -.d 
side  moments,  respectively.  If  the  roll  rate  is  very  small,  and  the  side 
moment  negligible,  then  the  roll  rate  dependent  terms  approach  zero  and  the 
damping  terms  must  also  approach  zero,  implying  that  (j)  is  near  zero.  For 
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the  purpose  of  obtaining  exact  solutions  for  it  is  useful  to  consider 
the  following  alternative  form  of  Equation  (2). 


where 


A  cos  0  +  B  sin  (Ji  =  C 


A  =  p 


C  =  p 


g  cos  Y 

J 

g  p  S  cos  Y  I 

r - ^ - J 

Sip  p  S  V  Cj^  p  S  V  g  Oj.  sin  Y  ] 

.  4  *  2in  V  ^  J 


C3„  p  S  d 


Solutions  to  Equation  (3)  exist  for 

<  A^  + 

For  flat  trajectories  (y  0)  the  latter  inequality  can  be  expressed  as 


"M  p  S  d^  V 


I,  P  S  V 


J- 


C5„  p  S  d 


g  p  S  d^ 


>7 


Examination  of  (4)  reveals  that  large  values  of  p,  or  Cg^  preclude  the 

~  p 

attainment  of  a  solution  for  while  large  damping  ,  makes  a  solution 

more  likely.  However,  if  is  small  and  can  be  negl|BCted,  there  may  be 

a  range  of  [  p  |  <  [  per  |  for  which  solutions  for  0  exist.  In  such 

cases,  we  will  define  as  the  critical  roll  rate.  At  roll  rates  less 
than  the  critical  rate, the  coning  rate  is  zero;  at  roll  rates  greater 
than  the  critical  roll  rate,  the  coning  rate  is  finite. 


Under  those  conditions  where  a  solution  for  <p  is  possible,  it  is  worth¬ 
while  to  examine  the  nature  of  the  solutions.  First,  consider  the  case 
whore  is  negligible  and  only  Cm.-,  is  important.  It  is  found  from 

Equation  (2)  that  the  variation  of  <|)  with  p  for  |  p  |  <  IPqi^!  tN^pically 

like  the  result  plotted  in  Figure  4,  which  was  computed  using  the  original 
estimate  for  in  Figure  3. 

P 


Figure  4.  Variation  of  Angle  of  Attack  Plane 

Orientation  with  Roll  Rate  for  S-Curve 
Bomblet 


A  significant  result,  as  cart  be  seen  fromJFigure  5,  is  that  as  the  roll 
rate  decreases  to  zero  the  solutiorts  for  <p  also  approach  zero.  In  the 
limiting  case  of  zero  roll  rate  the  bomblet  continually  flies  nose  up  with 
the  lift  vertical.  Further  analysis  of  Equation  (2)  shows  that  for  positive 
roll  and  positive  Magnus  the  nose  of  the  bomblet  points  to  the  left  of  the 
velocity  vector  (as  viewed  looking  in  the  direction  of  flight).  Similarly, 
the  nose  of  the  bomblet  points  to  the  right  when  the  magnus  moment  is 
negative. 
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v'.rtical 

plane 


Figure  5.  Description  of  Angle  of  Attack  Plane 

Orientation  and  Aerodynamic  Roll  Angle 


To  confirm  that  the  steady-state  solutions  are  indeed  stable  solutions 
(i.e.,  that  any  oscillation.  ,  is  damped)  a  perturbation  differential 
equation  is  vrritten  using  as  the  dependent  variable.  Analysis  of  the 
perturbaition  equation  shows  that  stable  solutions  exist  only  for 
-  Tr/2  <  <>  <  II/2.  Outside  this  region  the  A(|>  perturbations  are  unstable. 
This  implies  that  the  nose  of  the  bomblet  is  always  nose-up  with  respect 
to  the  trajectory  when  stable  zero-coning  solutions  exist. 

Now,  consider  the  case  where  the  classical  spin  dependent  magnus 
moment  is  negligible  and  the  side  moment  C_j^  is  present.  Equation  (4) 
implies  a  maximum  for  irrespective  of  the  trim  angle  of  attack. 

Consider  also  the  effect  of  when  C„  is  a  harmonic  function  of 
the  aerodynaiiiic  roll  angle,  $.  This  functioi^l  relationship  exists  for 
most  finned  bodies.  A  general  solution  in  the  form  of  Equation^(l)  is  not 
readily  attainable  because  of  the  complex  relaticnship  between  ^  and  $. 
However,  the  special  case  where  p  =  (Ji  =  0  is  worthy  of  examination,  since 
then  We  can  rewrite  the  contribution  of  ($)  as 
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dC 


SM 


d<I) 


A-j) 


when  the  bomblet  symmetry  planes  are  approximately  aligned  with  the  angle 
of  the  attack  plane. 


Examination  of  the  stability  criteria  for  perturbations  A(J)  about  (|>  = 
constant 'shows  that  for  stable  solutions 


d<I) 


(5) 


In  other  words,  the  angle  of  attack  plane  will  rotate  until  $  attains  a 
value  where  relationship  (5)  is  satisfied.  Six-degrees-of-freedom  motion 
computations  show  that  inequality  (S)  is  valid  even  when  the  spin  rate  is 
finite. 

Six-Degrees-of-Freedom  Motion  Simulations.  Six-degrees-of-freedom 
simulations  have  been  used  to  verify  the  zero  coning  theory,  to  examine 
the  initial  transient  motion,  and  to  assess  the  effect  of  the  highly  non¬ 
linear  aerodynamic  coefficients.  These  simulations  were  accomplished  for 
the  cylindrical  afterbody  S-Curve  bomblet  configurations  (Figure  1-a  and 
Table  I)  using  a  special  six-degrees -of  freedom  trajectory  program. 

The  dispenser  event  conditions  are  typical  of  subsonic  aircraft  weapons 
delivery. 


TABLE  I.  ASSUMED  PHYSICAL  CHARACTERISTICS  OF 
S-CURVE  AND  ROLL-THROUGH-ZERO  BOHBLETS 


Body  diameter  -  2  inches 

Body  length  -  8  inches 

Weight  -  1.51  pounds 

Ix  -  1.629  X  10-4  slug-ft^ 

I  -  1.06  X  10-3  slug-ft2 


iFor  six  fins,J  =  o,  Tr/6,  tt/S,  2  tt/S,  etc. 
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A  comparison  of  the  six-degrees-of-freedom  motion  with  that  predicted 
by  the  zero  coning  theory  is  illustrated  in  Figure  6.  The  coning  motion 
IS  represented  by  the  variable,  which  represents  the  orientation  of  the 
angle  of  attack  plane  with  respect  to  horizontal  (Figures).  Several 
arbitrary  initial  orientations  of  the  angle  of  attack  plane  are  assumed 
for  the  six-degrees-of-freedom  motion  calculations.  The  initial  angle  of 
attack  is  5  degrees.  For  these  calculations  only  the  spin  dependent  Magnus 
moment  is  assxjmed  to  be  present  Coi“iginal  estimate  curve  of  3),  and  the  fin 
cant  is  0.02  degree  .  This  cant  angle  results  in  a  steady-state  roll  rate 
less  than  the  predicted  critical  roll  rate. 


Initial  conditions 


V£  ■  600  ft/sec 


TIME  -  SEC 


Figure  6«  Motion  Historie.,  for  S-Curve  Bomblet  with 
0.02  Degree  Fin  Cant 
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As  can  be  seen  from  the  motion  history,  the  angle  of  attack  plane 
converges  to  a  fixed  orientation  after  about  six  seconds  of  flight, 
regardless  of  the  initial  orientation.  The  steady-state  value  of  0 
closely  approaches  the  theoretical  zero-coning  solution.  Note  that  the 
angle  of  attack  rapidly  approaches  the  nominal  trim  angle  of  attack  and 
that  the  transient  motion  is  heavily  damped. 

With  an  increase  in  the  cant  angle  to  0.1  degree,  the  angle  ^  does 
not  approach  a  steady-state  value,  but  continues  to  increase,  negatively. 
Figure  7  illustrates  this  motion,  which  can  be  compared  directly  with  the 
results  for  the  smaller  cant  angle.  Although  for  this  case  the  coning  rate 
is  finite,  nonetheless,  the  coning  rate  is  only  about  6  percent  of  the  roll 
rate. 


Initial  conditions 
Vg  “  600  ft/sec 


P 

RAO/SEC 


lOh 


typical  _ 
for  all  ^ 
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TIME  -  SEC 


Figure  7.  Motion  Histories  for  S-Curve  Bomblet 
with  0.1  Degree  Fin  Cant 


When  the  aerodynamic  coefficients  are  modified  to  include  roll  angle 
dependency  (and  fin-body  lateral  misalignment  and  static  unbalance. are 
also  introduced)  the  bomblet  dynamics  become  very  complicated.  For  example, 
when  the  moment  vector  representing  the  lateral  misalignment  is  aligned 
with  the  angle  of  attack  plane,  this  moment  acts  like  a  magnus  moment  and 
tends  to  cause  coning.  Interestingly,  this  effect  becomes  more  pronounced 
as  the  roll  rate  due  to  fin  cant  is  decreased. 
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Koll  tor<iucs  due  to  static  unbalance,  fin  cant,  and  the  effect  of 
aerodynamic  roll  angle  on  the  fins,  as  well  as  the  side  moment,  can  con¬ 
tribute  to  anomalous  rolling  motions  and  roll  lock-in  conditions.  Thirty- 
five  trajectories  and  motion  histories  were  computed  using  the  subsonic 
dispenser  event  conditions  given  to  Table  II  and  the  following  bomblot 
asymmetries: 

Fin  cant  =0.1  degree  (constant) 

Lateral  misalignment  =  0.3  degree  (standard  deviation) 

Static  unbalance  =  0.0016  inch  (standard  deviation) 


TABLE  II.  EVENT  CONDITIONS  FOR  REPRESENTATIVE 
DISPENSER  FUNCTION 

Subsonic  Delivery 

Ve  =  900  ft/sec 

Y  *45  deg 
h(HOB)  *  2000  ft 

Transonic  Delivery 

Ve  =1145  ft/sec 

Y  =  65.7  deg 
h(HOB)  *  2000  ft 

Supersonic  Delivery 

Ve  =  2000  ft/sec 

Y  *  50  deg 
h(HOB)  =  2000  ft 


Random  values  for  the  initial  angular  orientation  of  the  angle  of  attack 
plane,  fin  symmetry  planes,  static  unbalance,  lateral  misalignment,  and 
initial  pitch  disturbance  were  selected.  Out  of  the  thirty-five  trajectories, 
twenty-one  developed  a  roll  lock-in  type  motion  with  9  constant,  while  for 
the  remainder  4  varied  continuously.  For  those  that  locked-in,  tlie  lock-in 
angles  are  plotted  in  Figure  8,  along  with  the  variation  of  the  side  moment 
with  4.  As  can  be  seen,  lock-in  occurs  only  when  the  gradient  of  the  side 
moment  is  positive. 
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y///////  «9lons  of  roll  lock-ln 


Figure  8.  S-Curve  Bopblet  Roll  Lock-in  Characteristics 


With  zero  fin  cant  the  percentage  of  those  bomblets  experiencing  lock- 
in  increased  to  100  percent,  while  for  0.5  degree  cant  no  lock-ins  were 
encountered.  Further  note  should  be  made  of  the  fact  that  when  lock-in 
occurred,  the  angle  of  attack  remained  well  behaved,  and  for  the  configu¬ 
rations  considered,  there  was  no  noticeable  tendency  for  the  trim  angle  of 
attack  to  increase. 

While  the  flight  dynamics  of  the  S-Curve  type  bomblet  must  be  reason¬ 
ably  well  behaved,  the  essential  requirement  is  the  achievement  of  a  large 
and  reasonably  uniform  dispersion  pattern.  Thus,  the  integrated  effect  of 
the  motion,  as  related  to  the  deflection  of  the  bomblet  trajectories  is  the 
only  thing  that  really  matters.  Also,  even  though  large  dispersions  are 
desired,  some  bomblets  with  poor  dispersion  can  be  accepted  as  a  means  of 
filling  the  center  of  the  pattern.  This  leads  one  to  consider  the  dis¬ 
persion  from  a  purely  statistical  point  of  view.  The  Monte  Carlo  method 
is  well  suited  to  this  type  of  problem  and  its  use  for  dispenser  pre¬ 
diction  is  described  in  a  subsequent  section. 

Roll -Through-Zero  Bomblet 

Aerodynamic  Design  Considerations.  The  aerodynamic  design  of  a  roll- 
through- zero -type  bomblet  with  large  body-fixed  trimmed  lift  would  appear 
straightforward.  A  means  must  be  provided  for  trimming  the  bomblet  to  a 
large  angle  of  attack,  and  a  small  roll  torque  must  be  introduced  to 
reverse  the  initial  rolling  motion  of  the  bomblet.  Unfortunately,  the 
means  provided  for  trimming  the  bomblet  can  often  lead  to  sizeable  roll  and 
side  moments  which  are  dependent  upon  the  roll  orientation  of  the  bomblet. 
Thus,  under  transient  motion  conditions,  such  as  occur  immediately  after 
dispenser  opening,  these  induced  moments  can  completely  overpower  the  small 
roll  torque  introduced  for  reversing  the  roll. 

Six-degrees-of-freedom  motion  studies  of  a  high-lift  bomblet  con¬ 
figuration  with  small  wing-like  surfaces,  which  had  large  roll  dependent 
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moments,  showed  that  this  type  of  bomblet  has  very  erratic  roll  motion  and 
poor  dispersion.  Consequently,  an  alternate  bomblet  concept  similar  in 
external  geometry  to  the  S-curve  bomblet  (Figure  Ic)  was  selected  for 
investigation.  Two  of  the  fins  are  given  an  incidence  of  10  degrees  for 
the  purpose  of  trimming  the  bomblet  to  a  finite  angle  of  attack. 

Figure  9  illustrates  the  aerodynamic  characteristics  of  the  bomblet 
as  a  function  of  angle  of  attack  with  the  bomblet  in  the  trimmed  roll 
orientation.  The  variations  of  the  coefficients  with  roll  are  depicted 
in  Figure  10.  The  variations  of  the  forces  and  moments  with  roll  are 
sizeable  but  less  than  those  of  the  original  high-lift  configuration.  The 
loss  of  lift  due  to  the  aft  placement  of  the  aerodynamic  trim  surfaces  is 
apparent . 


a  -  OEG 


Cn. 


Figure  9.  Roll-Through  Zero  Bomblet  Aerodynamic 

Characteristics  as  a  Function  of  Angle  of 
Attack  (Trimmed  Roll  Orientation) 
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Figure  10.  Effect  of  Aerodynamic  Roll  Angle  on  Aerodynamic 
of  Roll -Through-Zero  Bomblet 

Particularly  noteworthy  from  Figure  10  is  the  fact  that  the  gradient 
of  the  roll  torque  with  respect  to  the  aerodynamic  roll  angle  at  the  trim 
orientation  is  positive  or  unstable.  Also  note  that  the  side  moment 
gradient  is  positive.  The  latter  has  a  very  significant  stabilizing 
influence  on  the  aerodynamic  roll  angle  by  nature  of  the  fact  that  the 
coning  motion  induced  by  tends  to  return  the  aerodynamic  roll  angle 

to  the  trim  orientation.  This  stabilizing  effect  would  become  de- 

stabilirir.g  if  the  aerodynamic  surfaces  used  for  trim  were  placed  forward 
of  the  center  of  gravity. 
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simplified  Motion  Theory.  A  simplified  closed-form  solution  is 
readily  obtains  for  the  motion  of  a  roll -through-zero-type  bomblct  during 
Its  dispersal  phase  of  flight.  If  induced  aerodynamic  effects  are  neglected 
and  we  assume  constant  velocity,  linear  aerodynamics,  a  constant  rol’  torque 
coefficient,  ,  and  trimmed  flight,  the  solution  of  equations-of-motion 
for  the  roll  o  and  lateral  translation  degrees-of-freedom  from  the  time  • 
at  which  p  =  0  to  time  ^  =  x  is 


Ve  +  '  *e'  -yf 

d  ""  2m’  ^ 


2ii  I. 


% 


C(x)  +  i  S(x) 


(6) 


where  C  (x)  and  S  (x)  are  Fresnel  integrals 

o 

X 


cos  ( — )  du 


S(x) 


=J  8in(  — 


-)  du 


and 


The  Fresnel  integrals  have  asymptotic  values  for 

<j)  -*■  00 

C  ix)  =  i 


such  that 


S  (X)  =  2 


C  (x)  +  iS  (x) 


(7) 
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The  same  result  as  equation  (7)  is  approximated  hy  (fr  =  ir/A,  where 

i  '  (8) 


C  {  +  iS(  i  ) 

V2  V2 


Therefore,  starting  from  zero  roll  rate,  the  lateral’ deflection  of  the  ’ 
trajectory  after  rolling  through  an  angle  $  =  tt/4  is  nearly  identical  to 
the  deflection  which  would  occur  after  a  large  number  of  roll  cycles.  ' 
If  roll -through-zero  occurs  and  there  exists  at  least  one-eighth,  of  a  roll 
cycle  before  and  after  the  roll-rate  reversal,  then  the  total  deflection 'of 
the  trajectory  is  approximated  by 


An  important  result  obtained  from  the  above  relations^iip  is  that  very  large 
trajectory  deflections  can  be  achieve^,  provided  that  a  small  roll  torque 
is  used  to  de-spin  the  bomblet.  ,  1  ■ 


An  estimate  of  the  dispersion  at  ground  level  can  easily  be  made,  if 
the  bomblet  has  a  large  ballistic  coefficient  such  that  the  flight  path  , 
following  roll -through-zero  is  relatively  straight.  The  along-range  and 
cross-range  dispersions  (R  and  r,  respectively)  due  to  Ay  are  simply 


!  (IQ) 


(11) 


R  = 


lA 

sin^Y 


AY 


r  = 


sin  Y 


AY 


Six-Degrees-of-Freedom  Motion  Simulations.  Figure  11  illustrates  ^ 
typical  set  of  motion  histories  for  the  roll -through -zero  bomblet,  based 
on  the  subsonic  event  conditions  of  Table  II.  The  data  are  selected 
from  one  of  a.  series  of  trajectories  comprising  a  Monte  Carlo  dispersion  . 
study.  This  particular  bomblet  has‘an  initial  negative  roll  irate  of  1.0 
rps,  a  mean  effective  fin  cant  of  0.115; degree,  and  a  eg  lateral  offset  of 
0.001  inch.^  ' 


^It  is  to  be  noted  that  for  this  particular  bomblet  the  maximum  induced 
roll  moment  due  to  eg  offset  is  about  '60  percent  of  that  due  to  the 
intentional  fin  cant.  '  1 
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Figure  11.  Motion  Histories  for  Roll -Through -Zero 
Bomblet 

Examination  of  Figure  11  shows  that  about  0.5  second  is  required  for 
the  bomblet  to  attain  zero  roll  rate  and  for  the  aerodynamic  roll  angle, 
to  stabilize  at  the  trim  orientation,  ♦  «  ir.  The  special  orientation 
of  the  angle  of  attack  plane,  angle  9  ,  remains  approximately  within  one 
quadrant  during  the  first  1.0  second  of  flight,  and  this  is  the  period 
during  which  the  trajectory  is  deflected.  The  trim  orientation  remains 
stable  at  =  TT  for  approximately  3.0  seconds  of  flight  (thus  substan¬ 
tiating  that  the  roll  torque  need  not  be  stable),  but  subsequently  the  spin 
rate  approaches  the  resonance  spin  rate  and  the  familiar  phase  shift  in 
the  trim  orientation  begins  to  occur. 

The  lateral  eg  offset  is  a  very  critical  factor  for  the  roll-through- 
zero-type  bomblet.  Because  the  aerodynamic  roll  angle  under  trim  conditions 
is  nearly  invariant,  the  induced  roll  torque  due  to  eg  offset  results  in 
random  variations  in  roll  torque  and  roll  rate.  If  these  become  large, 
the  roll:  rate  builds  up  to  values  which  approach  or  exceed  the  pitch 
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resonance  frequency.  Then  the  trim  phase  shift  causes  the  side  moment 
to  increase  with  possible  adverse  effect  on  the  dynamic  stability  such 
that  very  large  coning  angles  ensue. 

These  considerations  complicate  the  aerodynamic  design  since,  if  the 
designer  attempts  to  increase  the  pitch  frequency  to  provide  more  lateral 
eg  tolerance,  the  trim  moment  requirements  also  tend  to  increase.  This, 
in  turn,  leads  to  increased  induced  aerodynamic  moments  from  the  larger 
trim  surfaces. 


20 


SECTION  III 


MONTE  CARLO  DISPERSION  PREDICTION 


Procedure 


The  motion  of  an  actual  S-Curve  bomblet  or  roll-through- zero  bomblet 
differs  considerably  from  that  of  an  ideal  bomblet  due  to  configurational 
asymmetries  and  departures  from  complete  rotational  symmetry.  In  addition, 
the  initial  motion  imparted  to  the  bomblets  at  dispenser  opening  and  during 
the  subsequent  period  of  cluster  breakup  tends  to  be  stochastic  in  nature. 

The  Monte  Carlo  method  provides  a  means  for  ascertaining  the  dispersion 
characteristics  of  a  bomblet  cluster  by  computation  of  the  probable  motion 
and  impact  points  of  a  large  number  of  bomblets.  This  was  accomplished  by 
the  use  of  statistical  values  for  the  initial  motion  parameters  and  bomblet 
configurational  asymmetries.  For  each  simulated  dispenser  event,  from  30 
to  as  many  as  100  six-degree-of-freedora  trajectories  are  computed  to 
determine  the  probable  dispersion  and  impact  pattern  of  the  bomblet  cluster. 

The  Monte  Carlo  simulations  described  in  this  report  have  been  com¬ 
puted  primarily  for  the  representative  subsonic  dispenser  event  conditions 
described  in  Table  II.  The  bomblet  physical  parameters  correspond  to  the 
data  of  Table  I.  The  bomblet  size  and  inertial  properties  are  repres,entative 
of  a  hypothetical  configuration. 

Initial  Motion  Parameters  ^  , 

The  initial  motion  parameters  are  evaluated  in  accordance  with  the 
probable  time  of  bieak-out  from  the  cluster,  defined  as  tb^.  At  this  instant 
the  bomblet  is  no  longer  influenced  by  the  multi -component  flow  effects 
which  exist  because  of  bomb let -bomb let  and  bomblet -dispenser  interference 
effects.  The  statistical  values  of  the  motion  parameters  at  break-out  are 
determined  empirically  from  flight  test  data  or  suitable  assumptions. 

The  probable  cluster  break-out  time  is  established  using  a  special 
distribution  function.  Figure  12  shows  a  typical  curve  of  cumulative 
frequency  versus  cluster  break-out  time  based  on  experimental  data. 

The  bomblet  velocity  at  ti,  is  assumed  to  be  comprised  of  two  components 
a  component  in  the  mean  °  direction  of  motion  of  the  cluster,  and  a 

radial  velocity  component,  with  a  Gaussian  magnitude  (but  not  necessarily 
zero  mean)  and  uniformly  random  orientation.  The  longitudinal  component  of 
velocity  depends  upon  the  event  velocity  and  the  cluster  deceleration.  The 
bomblet  cross-angular  velocity  is  assumed  to  be  Gaussian  in  magnitude  and 
uniformly  random  in  direction.  The  initial  angle  of  attack  is  arbitrarily 
selected  as  zero  because  the  cross-angular  velocity  establishes  the  magni¬ 
tude  of  the  first  oscillation  peak. 
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ACCUMULATIVE 

PROBABILITY 


Figure  12.  Statistics  of  Bomblet  Cluster  Breakup 


For  the  subsonic  event  conditions  given  in  Table  II,  the  standard 
deviations  for  the  initial  radial  velocity  and  initial  cross-angular 
velocity  are 


Vj^  a  20  ft/sec 
q*  =  35  rad/sec 
Configurational  Asymmetries 

The  physical  variations  in  bomblet  configuration  considered  in  the 
present  Monte  Carlo  simulations  are  fin  cant,  longitudinal  misalignment, 
and  lateral  eg  offset.  The  latter  two  parameters  are  considered  to  be 
Gaussian  in  magnitude  and  uniformly  random  in  orientation.  In  general,  the 
fin  cant  has  a  non-zero  mean  and  Gaussian  distribution.  Since  there  are  no 
S-Curve  or  roll -through-zero  bomblets  in  production,  estimates  for  the 
configurational  asymmetries  are  based  on  other  representative  ordnance. 

Values  for  the  standard  deviations  of  fin  cant  error,  longitudinal 
misalignment  and  lateral  eg  .offset  are  selected  on  the  basis  of  fabri¬ 
cation  accuracies  which  are  representative  of  small  ordnance.  Of  these, 
the  standard  deviation  of  eg  lateral  offset  is  most  sensitive  to  the  par¬ 
ticular  fabrication  and  loading  techniques  ured. 

For  most  of  the  Monte  Carlo  simulations  the  standard  deviation  of 
eg  lateral  offset  is  0.0016  inch,  although  a  value  as  large  as  0.0400  inch 
is  considered  for  the  S-Curve  bomblet. 
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The  standard  deviation  value  for  effective  fin  cant  error  is  dependent 
upon  the  number  of  fins  since,  if  the  variance  of  all  the  fins  is  identical, 
the  standard  deviation  of  roll  torque  is 

standard  deviation  C.=  — — — 

*  -\/^ 

where  Op  is  the  standard  deviation  of  error  of  a  single  fin.  For  a  six- 
fin  configuration,  the  standard  deviation  of  effective  cant,  5g££.  is 

=  0*41  Op.  The  effective  longitudinal  misalignment  of  the  fin 

assembly  is  similarly  defined  by  S  *  ^eff  For  six-fins  and  neglecting 
fin-fin  interference,  the  standard  deviation  of  for  longitudinal  mis¬ 
alignment  is  =  1.12  Op.  This  is  combined  with  the  body  misalignment 

on  an  RSS  basis.  In  the  present  studies  a  value  of  0.2  degree  is  used  for 
both  Op  and  the  standard  deviation  of  body  misalignment,  based  on  an 

averaging  of  several  data  sources. 
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SECTION  IV 


MONTE  CARLO  IMPACT-DISPERSION  PATTERNS 

S-Curve  Bomblet 

A  computed  Monte  Carlo  dispersion  pattern  for  the  type  a  S-Curvo 
bomblet  configuration  is  depicted  in  Figure  13,  For  this  simulation,  the 
bomblets  have  a  mean  fin  cant  of  zero  degrees,  and  a  standard  deviation 
of  eg  lateral  offset  of  0.0016  inch. 


i-ROSS-RAKCE  -  FT 


Hole:  subsonic  dispenser 
event  conditions 
(Table  I) 


Figure  13.  Monte  Carlo  Impact-Dispersion  Pattern  for  S-Curve  Bomblet 

Statistical  analysis  of  the  impact  pattern  shows  that  the  standard 
deviations  of  along-range  and  cross-range  dispersion,  with  respect  to  the 
mean  center-of-impact,  are  ±  246  feet  and  ±  185  feet,  respectively. 


24 


The  maximum  dispersion  is  about  2.3  times  greater  than  the  standard 
deviation  values,  which  emphasizes  the  fact  that  the  impact  points  are 
concentrated  toward  the  center  of  the  pattern. 

Accuracy.  To  obtain  an  approximate  indication  of  the  impact  pattern 
statistical  repeatability,  a  set  of  33  trajectories  was  rerun  several  times 
using  the  same  statistical  distribution  functions  but  with  different  random 
numbers.  It  was  determined  that  for  a  Monte  Carlo  sample  size  corresponding 
to  33  trajectories,  the  standard  deviations  of  along-range  and  cross-range 
dispersion  are  accurate  to  within  about  ±25  feet  at  the  90  percent 
confidence  level . 

Comparison  with  Theory.  A  theoretical  dispersion  boundary  is  plotted 
in  Figure  13  for  purposes  of  comparison  with  the  Monte  Carlo  results.  The 
theoretical  boundary  represents  the  dispersion  which  would  be  achieved  with 
zero  coning.  In  determining  the  zero  coning  boundary,  an  average  time  for 
cluster -breakout  was  assumed.  As  can  be  seen,  several  impact  points  approach 
the  zero-coning  boundary,  and  this  number  would  be  expected  to  increase  in 
proportion  to  the  total  number  of  trajectories  simulated.  Thus,  the  zero 
coning  theory  provides  a  close  approximation  of  maximum  pattern  size. 

Effect  of  Static  Unbalance.  Monte  Carlo  impact  patterns  have  been  com- 
puted  for  various  arbitrary  values  of  the  standard  deviation  of  lateral  eg 
offset.  The  trend  of  the  maximum  dispersion  and  the  standard  deviation  of 
dispersion  with  the  standard  deviation  of  lateral  eg  offset  is  depicted  in 
Figure  14,  There  is  a  decrease  in  both  the  maximum  and  standard  deviation 
values  of  dispersion  with  increasing  eg  offset.  The  dat^  show  that  a  eg 
offset  tolerance  of  0.01  inch  can  be  applied,  if  a  10  percent  reduction  in 
cross-range  dispersion  can  be  accepted.  It  should  be  noted  that  the 
sensitivity  to  static  unbalance  increases  appreciably  when  a  larger  trim 
angle  of  attack  is  selected. 

Effect  of' Fin  Cant.  It  has  been  theorized  that  the  coning  motion  due 
to  the  body-fixed  roll -dependent  side  moments  would  be  less  if  an  inten¬ 
tional  roll  rate  was  introduced  such  that  the  effect  of  these  moments  on 
the  flight  path  would  tend  to  be  averaged.  Dispersion  patterns  were  com¬ 
puted  with  a  mean  fin  cant  of  O.S  degree  for  a  range  of  static  unbalance. 

The  results  are  compared  with  the  zero  cant  data  in  Figure  15.  With  the 
intentional  cant  angle,  there  is  a  large  increase  in  the  standard  deviation 
values  of  along-range  and  cross-range  dispersion,  thus  proving  the  effective¬ 
ness  of  this  bomblet  modification. 

Effect  of  Initial  Motion  Parameters.  The  sensitivity  of  the  S-Curve 
bomblets  to  the  magnitude  of  the  initial  motion  disturbance  has  been  investi¬ 
gated  using  the  Monte  Cavlo  method.  Doubling  the  initial  tumble  rate  to 
70  rad/sec  had  negligible  effect  on  the  dispersion,  even  though  the 
magnitude  of  the  first  oscillation  peak  was  as  large  as  90  degrees. 
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Figure  14.  Effect  of  Static  Mass  Unbalance  on  S-Curve 
Bomblet  Dispersion  (Monte  Carlo  Results) 

““  W«n  fin  c»nt  •  0.5  deo 
Clean  fin  cant  •  0 

DISPERSION  -  FT 
400 


300 


200 


1JJ 


U 

0.0001  0.001  0.01  0.1 
SJANOARD  DEVIATION  OF  C.G.  OFFSET  -  INCHES 

Figure  15.  Effect  of  Mean  Fin  Cant  Angle  on  S-Curve 
Bomblet  Dispersion  (Monte  Carlo  Results) 
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Effect  of  Increasing  Trim  Angle  of  Attack.  Increasing  the  bomblet 
nominal  trim  angle  of  attack  from  10  degrees  to  16  degrees  resulted  in  a 
dispersion  increase  (standard  deviation)  nearly  proportional  to  the  ratio 
of  trim  lift  coefficient.  These  results  were  achieved  in  spite  of  the 
increased  nonlinearity  of  the  aerodynamic  coefficients  at  the  larger  angle 
of  attack. 

Transonic  and  Supersonic  Event  Conditions.  For  transonic  delivery 
conditions,  calculations  have  shown  the  boattail  afterbody  configuration 
(Figure  lb)  to  have  more  desirable  dispersion  characteristics  than  the 
cylindrical  afterbody  bomblet,  even  though  the  trim  angle  of  attack  is 
less  than  that  of  the  cylindrical  afterbody  for  Mach  numbers  between  1.0  and 
1.5.  Monte  Carlo  dispersion  data  for  the  boattail  afterbody  bomblet  are 
summarized  in  Table  III  for  the  representative  transonic  and  supersonic 
event  conditions  as  given  in  Table  II. 


TABLE  III.  TRANSONIC  AND  SUPERSONIC  DISPERSION 
CHARACTERISTICS  OF  S-CURVE  BOMBLET 
WITH  BOATTAIL  AFTERBODY 


Standard 

Standard 

Deviation 

Deviation 

Range 

Cross-Range 

Dispersion 

Dispersion 

(Feet) 

(Feet) 

Transonic  Event 

202 

183 

Supersonic  Event 

959 

470 

For  transonic  delivery,  the  impact  dispersion  is  comparable  to  that 
for  the  cylindrical  afterbody  bomblet  delivered  at  high  subsonic  velocity. 

In  contrast,  for  supersonic  delivery,  the  range  dispersion  is  increased  by 
a  factor  greater  than  four  and  the  cross-range  by  a  factor  greater  than  two. 

Roll-Through-Zero  Bomblet 

A  Monte  Carlo  dispersion  pattern  for  the  roll -through -zero  type  is 
illustrated  in  Figure  16.  The  standard  deviation  values  of  range  and  cross¬ 
range  dispersion  are  207  feet  and  159  feet,  respectively.  The  dispersion 
pattern  has  a  maximum  diameter  of  about  650  feet.  For  this  simulation  the 
standard  deviation  of  bomblet  static  unbalance  is  0.0016  inch. 
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Figure  16,  Monte  Carlo  Impact  Dispersion  Pattern  for  Roll-Through- 
Zero  Bomb let 


Comparison  with  Theory.  Figure  16  also  depicts  the  theoretical  dis¬ 
persion  boundary  obtained  using  equations  (9),  (10),  and  (11)  with  the  roll 
torque  coefficient  in  equation  (9)  corresponding  to  the  mean  fin  cant  angle. 
There  is  excellent  agreement  between  the  theoretical  dispersion  boundary  and 
the  standard  deviations  of  range  and  cross-range  dispersion  computed  from 
the  Monte  Carlo  impact  pattern.  However,  using  the  mean  fin  cant  angle, 
the  maximum  dispersion  is  significantly  underpredicted,  as  would  be  expected. 

Comparison  with  S-Curve.  A  comparison  of  the  impact  pattern  data  for 
the  S-Curve  "and  roll -through- zero  bomblets  (standard  subsonic  event  conditions) 
shows  that  the  ratios  of  the  standard  deviation  values  of  dispersion  for  the 
two  types  of  bomblets  are  very  nearly  proportional  to  the  respective  trimmed 
lift  coefficients.  This  comparison  is  based  onanC-Curve  bomblet  config¬ 
uration  with  zero  mean  fin  cant.  For  supersonic  event  conditions,  the 
performance  of  the  roll -through-zero  bomblet  is  poor  compared  to  the  S-Curve 
bomblet,  because  the  roll-through-zero  bomblet  dispersion  does  not  benefit 
from,  increased  dynamic  pressure,^ 

®  Equation  (9)  shows  that,  theoretically,  the  roll-through-zero  bomblet 
dispersion  is  independent  of  velocity. 
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SECTION  V 


COMPARISON  OF  S-CURVE  AND  MAGNUS  ROTOR  BCWBLETS 


One  of  the  advantages  of  the  S-Curve-type  bomblets  is  that  the  aero¬ 
dynamic  dispersal  force  is  generated  immediately  upon  bomblet  ejection 
into  the  air  stream.  Further,  the  dynamic  pressure  is  large  during  the 
initial  phase  of  flight  and  this  aids  in  the  rapid  dispersal  of  the  bomblet 
cluster.  In  contrast,  Magnus  rotor  bomblets  must  undergo  a  period  of  spin- 
up  during  which  the  aerodynamic  lift  is  small.  The  lift  will  remain  small 
until  the  magnus  rotor  decelerates  to  below  approximately  0.5  Mach  number. 

Figure  17  compares  typical  S-Curve  and  magnus  rotor  trajectories 
for  the  representative  subsonic  dispenser  event  conditions  of  Table  II. 

The  more  rapid  dispersion  of  the  S-Curve-type  bomblet  is  clearly  evident. 
The  better  performance  of  the  S-Curve-type  bomblet  can  also  be  translated 
into  a  lower  height  of  burst  (HOB).  For  example,  assuming  the  same  impact 
dispersion,  the  S-Curve  bomblets  can  be  deployed  at  a  650-foot  lower 
altitude. 
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Comparison  of  S-Curve  Bomblet  and  Magnus  Rotor 
Tra j  ectories 
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SECTION  VI 


CONCLUSIONS 


On  the  basis  of  analytical  studies,  including  Monte  Carlo  impact 
dispersion  predictions,  both  the  S-Curve  and  roll -through-zero-type  bomblets 
are  found  to  have  good  potential  for  cluster  dispenser  and  cluster  warhead 
applications.  These  bomblets  must  be  very  carefully  designed  to  insure 
that  the  proper  aerodynamic  characteristics  are  obtained  and  attention 
must  be  given  to  the  allowable  configurational  asymmetries. 
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